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Abstract 

A model is presented that fits the quark and lepton masses and 
mixings wherein five dimensionless parameters and a phase account 
for fifteen dimensionless observables. Among these are the Wolfenstein 
parameters p and 77, the v e — z/„ and v e — v T mixing angles which are 
predicted to be small and comparable while the — v T mixing angle is 
predicted to be large. The model is based on supersymmetric 5O(10) 
with the form of the mass matrices motivated by simplicity at the 
level of grand unification. 



PACS numbers: 12.15.Ff, 12.10.Dm, 12.60.Jv, 14.60.Pq 



E-mail: albright@fnal.gov; 2 E-mail: smbarr@bartol.udel.edu 



In an earlier paper with Babu M, the authors have proposed a model 
based on supersymmetric SO (10) for the masses and mixings of the heavier 
two families of quarks and leptons. One of the guiding principles in con- 
structing the model was that the mass terms should be of a kind that can 
arise simply in SO (10) with small representations of Higgs fields 0. Two 
critical points emerged from that work: namely, that it is quite natural in 
unified theories for there to be flavor-asymmetric contributions to the fermion 
mass matrices, and that these can explain some of the puzzling features of 
the spectrum. In particular, it was shown that a flavor-asymmetric term can 
account in a simple way for the apparent fact that the — v T mixing is very 
large compared to V c b- 

The present work is an extension of that earlier model to include the first 
family. We find that it is only necessary to introduce two new parameters, 
one of which has a complex phase, in order to account for seven quantities 
that involve the first family. In doing so, several interesting predictions result, 
including the values of the Vz — v^ and v e — v T mixing angles, and a relationship 
between the real and imaginary parts of V u b, i.e., between the Wolfenstein 
parameters p and rj. Altogether, as will be seen, the model leads to nine 
predictions. 



The mass matrices. 

The mass matrices proposed have the following forms: 
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where U, D, L, and N stand, respectively, for the up-quark, down-quark, 
charged-lepton, and the neutrino Dirac mass matrices. The zero superscripts 
stand here and throughout for quantities defined at the grand unification 
(GUT) scale. The convention we use is that the left-handed fermions multiply 
these matrices from the right, while the right-handed fermions multiply them 
from the left. (Thus, it is the entries below the main diagonal that control 
the quark and lepton mixings that are observed.) 
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It is immediately apparent that the matrix U° is rank-2, corresponding 
to a massless u quark. We imagine that some small higher-order effects 
give u a mass of several MeV. The fact that m u must arise at higher order 
than ma and m e is not a flaw but a qualitative success of the model. The 
point is that a several-MeV w-quark mass corresponds to m°/m° « 10~ 5 . 
This is very much smaller than the corresponding ratios mj/m" « 10~ 3 and 
m° e /m® « 0.3 x 10~ 3 , and so it is reasonable to suppose that it is a higher- 
order effect. 

By field redefinitions one can make all the elements of the matrices in Eq. 
(1) real except for 5'. (There can also be a relative phase initially between 
a and e, which was called a in [0. After phase redefinitions, however, this 
parameter survives only in the entries D® 3 and Lg 2 , which can be approxi- 
mated by a + | cos a and a + ecosa, respectively, since as will be seen later 
e/cr is a small parameter. As cos a only appears in terms subleading in e/cr, 
it has only a few-percent effect on predictions. In |I| it was found that the 
fits prefer a = 0, so we will assume that value henceforth in this paper.) Ef- 
fectively, then, the mass matrices in Eq. (1) contain five real dimensionless 
parameters and one phase: (Mu/Md), &, e, 8, \5'\, and arg<5'. 

We will first describe the kinds of 5*0(10) operators and diagrams that 
give rise to the entries shown in (1), and explain the rationale for the those 
entries. Then we shall proceed to discuss the fit to the data and the predic- 
tions of the model. 

The origin and rationale of the matrices. 

The (33) entries of the mass matrices contribute to the masses of the 
heaviest family of quarks and leptons. As is true of most grand unified 
models of fermion masses, these entries can come from the minimal kind of 
Yukawa term. Since ours is an SO (10) model, these entries arise simply from 
the operator (163163)10^. Here and throughout, the subscript L H' stands 
for a Higgs multiplet, and a numerical subscript stands for the family number 
of a quark/lepton multiplet. The parameters Mjj and Mb in Eq. (1) come 
from the SU(2)l x [/(l)y-breaking vacuum expectation values (VEVs) in 
5(10//) and 5(10#), respectively. We use (p(q) to represent a p multiplet 
of SU(5) contained in a q multiplet of SO (10). 

While the masses of the third family satisfy the minimal SU(5) (and 
5*0(10)) relation = mf, the middle family does not satisfy the corre- 
sponding m' = m[J relation. As is well known, however, a good fit is given 
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by the Georgi-Jarlskog relation J|, m° s = m°/3. Perhaps the simplest way 
to get a factor of 1/3 into the elements of the mass matrices is through 
an adjoint Higgs multiplet, 45#. Such a Higgs multiplet is needed in any 
case for the breaking of SO (10) down to the standard model group. In the 
simplest breaking scheme that naturally solves the "doublet-triplet splitting 
problem", the so-called Dimopoulos-Wilczek or "missing VEV" mechanism 
[|j], there is an adjoint Higgs multiplet whose VEV is proportional to the 
SO (10) generator B — L. This is in fact just what is needed to get a relative 
factor of 1/3 between entries in the quark and lepton mass matrices. 

The lowest- dimension Yukawa operators that contains the adjoint Higgs 
multiplet are products of the form 16jl6jlO#45H. There are two ways to 
contract the SO (10) indices, but it is easy to show that if (45#) oc B — L, 
then only the flavor- antisymmetric piece contributes to the mass matrices. 
This accounts for the form of the "e" entries in Eq. (1). Note that these 
appear with a factor of 1/3 in the quark matrices and 1 in the lepton matrices, 
reflecting the generator B — L. One also sees that because the operator being 
discussed contains the same 10# that appears in the Yukawa operator for 
the (33) entries, the e entries of U° and N° are Mu/Md times those of D° 
and L°. 

The operator I62I63IO//45H needed for the e entries arises in a sim- 
ple fashion by integrating out a superheavy (i.e. GUT-scale) family plus 
conjugate-family pair: 16 + 16 and is naturally of order 1/M G as required. 
This is shown in Fig. 1. Such an operator could also arise from Planck-scale 
physics, but it would then be too small. 

The next pieces of the mass matrices to be discussed are the "cr" entries. 
These have several rationales, and are in fact the key to the whole model. We 
shall first explain how such entries arise, and then why they are important 
for getting a realistic model. 

The feature of the a entries that distinguishes them from what is found in 
most published models is that they are neither flavor-symmetric nor flavor- 
antisymmetric, but rather "lopsided": a appears in D° 23 but not in £>° 2 . 
Though it is not common to consider such lopsided forms, there is really no 
reason from the point of view of grand unification not to do so. In fact, it 
may well be the neglect of such possibilities that has handicapped the search 
for realistic quark and lepton mass matrices. 

The a entries arise very simply from integrating out an SO (10) vector 
10 of quarks and leptons as shown in Fig. 2. The resulting operator is of 
the form [162l6#] 10 [16316^] io- The notation [...]io indicates that what is 
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contained in the brackets is contracted into a 10 of SO (10) as a result of the 
virtual 10 fermions appearing in Fig. 2. The Higgs multiplet denoted 16h 
is assumed to develop a GUT-scale VEV in the 1(16) direction, while the 
other spinor Higgs, denoted 16' H , gets a weak-scale SU (2) ^-doublet VEV in 
the 5(16) direction. Thus, in SU(5) notation, the operator from Fig. 2 gives 
(5 2 10 3 )(5^)(1 H ). 

One observes that since the left-handed down-quarks are in 10's of SU(5), 
not 5's, this operator contributes to D23 and not D^- On the other hand, 
since the left-handed leptons are in 5's, not 10's, the operator contributes 
to L 32 and not L23. As we shall see, this fact that the a entry is transposed 
in L relative to that in D is very important for understanding the large 
Vjx — v T mixing that may have been seen in the atmospheric neutrino data. 
This transposition is a consequence of the group theory of SU (5) which also 
implies that D = L T in the minimal SU(5) model. 

The upshot is that integrating out 10's of fermions in 5*0(10) can quite 
typically give just such lopsided entries as the a entries in Eq. (1). Note that 
the a entries only appear in D and L, and not in U and N, following from 
the fact that 16'^ contains only 5 and not 5 of SU(5). 

Now, why should such lopsided entries help to get a realistic model? It 
turns out that these a entries play at least six vital roles: 

(i) In the first place, it is easy to see from the seesaw-type formulas, that 
setting a = in Eq. (1) would give m"/m" = (e/3) 2 and m^/mj = e 2 . Thus 
one would get a ratio of ratios of 1/9 instead of the desired Georgi-Jarlskog 
value of 1/3. If, however, a is present and is large compared to e, then one 
has m°/m° ~ cre/3, and m"/mj ~ ere, so that the desired Georgi-Jarlskog 
factor emerges. 

(ii) Without the er entries, the (23) blocks of D and U would be propor- 
tional, and consequently there would be no KM mixing for those families, 
that is V c b would vanish. 

(iii) Moreover, the lopsided er entries explain a peculiar fact about the 
magnitude of If _D 2 3 and D 32 were of equal magnitude, as is the case 
with the purely flavor-symmetric or antisymmetric forms usually assumed, 
then V c b would get a contribution from D that is ^/m s /m{,, which is much too 
big (by about a factor of 4). By having the large lopsided contribution with 
a > e, one finds D 32 <^ D23, which leads to V c b ^Jm s /mb, as observed. 
This cures one of the main problems with textures of the Fritzsch type || 
as noted previously by Branco and Silva-Marcos in a phenomenological 
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study of the quark sector. Other authors |7[] have also recently appreciated 
this fact in the context of SU(5) models. 

(iv) Because a naturally only appears in D and L but not in U, one 
explains why m c /m t <C m s jm^ and m^/m T , as actually observed. In the 
absence of the a entries, one would obtain the minimal 50(10) prediction 
m°/m° = m®/m®, which is wrong by about a factor of 20. 

(v) The fact that the large a entry appears transposed in L relative 
to D means that it makes the (32) mixing of the left-handed leptons large 
compared to the corresponding mixing of the left-handed quarks. One thus 
explains why — v T mixing is much larger than V c b, which is one of the great 
puzzles in understanding the atmospheric neutrino oscillations. 

(vi) Finally, we shall see later that the a entries explain why the v e — 

mixing angle is significantly smaller than ^rrie/m^, assuming the small-angle 
MSW solution of the solar neutrino problem. 

What remains to be explained about the mass matrices in Eq. (1) are the 
entries 8 and 8'. These arise very simply from the diagrams in Fig. 3. Note 
that because the effective operators are again of the form 16jl6jl6'^16#, 
as for the a entries, these diagrams only contribute to D and L, and not 
to U and N. Note also that because they come from integrating out a 10 
of Higgs, they are flavor-symmetric. These contributions are assumed to be 
much smaller than the others. 

The fact that the above terms appear and no others presupposes that 
some Abelian symmetries exist that distinguish among families. Such sym- 
metries have been used by many authors. For an interesting recent publica- 
tion with a number of references cited, see ||. We note that realistic models 
involving such symmetries can be constructed giving rise to the matrices in 
Eq. (2). Much of the structure of such a model has already been described 
in [l] and 0. In addition to what was presented there, there have to be 
50(10) 10's of Higgs fields to be integrated out to give the 8 and 8' entries, 
as in Fig. 3. Such fields and the necessary couplings can be introduced, with 
appropriate Abelian symmetry, without destabilizing the gauge hierarchy or 
creating other difficulties 0. 

The predictions of the matrices. 

In expressing the predictions of the model, it proves convenient to use in- 
stead of the parameters 8 and 8', the linear combinations t R = 8y/a 2 + 1/ (cre/3) 
and ti = (8 — a 8') /(cre/3). Essentially, t^ is the Cabibbo angle, and tR is 
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the corresponding mixing angle of the right-handed d and s quarks. That 
is to say, if one diagonalizes the (23) block of D° in Eq. (1), and calls the 
resulting matrix D ' , then t L = D%[/D% 2 and t R = D^/D^. In the phase 
convention being adopted, t R is real and has a phase that shall be called 
9. 

The parameters, then, are Mjj/Me>, a, e, t R , and 0. In terms of these 
one has the following relations, which are given only to leading order in the 
small parameters e, t R , and t^: 
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Predictions involving only the second and third families were given already 
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in [0, where a fit was done using more exact expressions, and taking into 
account renormalization-group running from the unification scale down to 
low scales as well as radiative corrections to the heavy quark masses coming 
from supersymmetry. From m^/m T and V c b which have small experimental 
uncertainties, e and a, which was previously called p in [0, were extracted; 
cf., Eqs. (2e) and (2i). It was found that the best fit was given by 

e S 0.139, 

(3) 

a S 1.81. 

From Eq. (2c) one can then obtain a prediction for the charm quark mass. 
In [0 it was found that m c {m c ) = 1.10 ± 0.1 GeV, in satisfying agreement 
with the experimental value 1.27 ± 0.1 GeV. 

From Eq. (21) one sees that the model predicts a large, and in fact a 
nearly maximal, mixing of and v T . The element {U^l, which describes 
this mixing, receives a large contribution from diagonalizing the mass matrix 
of the charged leptons, L. This large contribution is well determined in the 
model, since L is known, and comes out to be approximately a j\/a 2 + 1 = 
0.87. However, [U^l also receives a small contribution from diagonalizing the 
neutrino mass matrix, which has the usual seesaw form M v = — jV^M^iV 
[[UJ . While N is known in the model, Mr is not — at least not without some 
further assumptions. Nevertheless, one can say that the contribution to \U^\ 
from diagonalizing M v is of order e, following from the form of N. Though 
relatively small, this contribution makes impossible a precise prediction of 
the v lx —p T mixing angle, though it is possible to say that it is nearly maximal. 
This source of uncertainty does not affect the predictions for the other two 
lepton mixing angles, however, because of the fact that N, and therefore 
M„, have zeros in the first row and column. These zeros are doubtless lifted 
by the higher order effects that generate m u , but it is easy to see that this 
should have a negligible effect on {7 e2 and U e z. 

The other predictions for the heavy families are that m° = m° T from Eq. 
(2b) and m°/m° = |m°/m° from Eqs. (2d) and (2e). Both of these are 
well-known relations, and work well in the context of supersymmetric grand 
unification. The foregoing four predictions were thoroughly discussed in [0 . 

Let us turn now to the eight quantities that involve the first family: m e , 
rrid, m u , U e2 , U e3 , V us , Re(V u b), and Im(\4t,). These depend on the three 
additional parameters tn, \ti\, and 9. The approximate magnitudes of tu 
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and | t^l can be estimated in a simple way, as follows. Comparing the well- 
known empirical relation sin^c = ^m^Jm s with Eqs. (2j) and (2g), one sees 
immediately that ~ tR ~ sin^. This enables one to read off an estimate 
of the Wolfenstein parameters [JTTJ] p and r] from Eq. 2. To leading order in 
the Wolfenstein parameter A, one has that p + ii] = (y u }jV cs /V c iy us )* . But 
the value of that invariant is given in this model (to leading order in e, tR, 
and ti) by 
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See Eqs. (2i), (2j), and (2k). Since ~ tR and o ~ V3 one has, roughly, 
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corresponding to a circle in the p — i] plane. 



To get the predictions more exactly, it is best to use V us and m e /m M , which 
have the smallest experimental uncertainties, to determine tR and \tt\] cf., 
Eqs. (2h) and (2j). Employing more exact expressions than given in Eqs. 
(2), one finds that 



t R 



0.239 
0.185. 



(5) 



The phase of ti, denoted 6, is so far undetermined. Substituting these values 
into Eq. (4), and taking into account the corrections to that equation that 
are higher order in e, tR, and t^, one finds 



p + ir] = 0A7e l ° - 0.29. (6) 

This circle, plotted as the bold dashed arc in Fig. 4, slices nicely through 
the presently allowed region . 

The prediction for the v e — mixing, cf. Eq. 2(m), is that 

sin 2 26 efl = 4 \U e2 \ 2 (1 - \U e2 \ 2 ) = 6.6 x lO" 3 . (7) 

This is well within the present limits 3 x 10 -3 < sm 2 29 etl < 8 x 10~ 3 for 
the small angle MSW solution to solar neutrino problem. The crucial role of 
the parameter a should be noted. In its absence one would find the familiar 
type of prediction, U e 2 — ^Jm e /m^, which is much too large. However, since 
there is a large intermingling of the second and third families caused by the 



9 



parameter a — the intermingling that causes the large — v T mixing — the 
angles \J e 2 and \J e ^ get intermingled as well. This leads to the suppression 
of U e 2 by a factor of cos6^ T ~ 0.7, which brings it into the experimentally 



allowed range 12 



This same intermingling gives rise to a distinctive prediction of this model, 
namely that 



'e3 



tan0 MT |[/ e2 | . (8) 



Using more exact expressions, one finds that 

sin 2 29 eT ^ 1.1 tan 2 6^ sin 2 29 e ^. (9) 

The remaining two predictions of the model are that m u ~ 0, Eq. (2f), 
which has been discussed already, and that m° d /wP s = Qm^/m^o, Eqs. (2g) 
and (2h), the other Georgi-Jarlskog relation, which gives m^/m s = 1/23 in 
excellent agreement with the current- algebra result. 

A comment is in order about the robustness of these predictions. In ob- 
taining the value of £r, we used m e /m^ on the grounds that the experimen- 
tal uncertainty in this quantity is negligible. However, one might question 
whether Eq. (2h) is trustworthy. The point is that since some higher order 
effects are presumed to generate a mass for the u quark, the same effects 
are likely to make a small contribution to the masses of d and e as well, 
though that is not necessarily the case. For example, m u could arise from 
~ 10~ 5 M[/, and a comparable entry in L of order L\ Y ~ 10~ 5 Md would 
affect m e at the 3% level. This would change £r, and thus U e 2, by the same 
amount. The effect of such uncertainties on the prediction for the Wolfen- 
stein parameters (p, rj) is relatively small. The prediction in Eq. (9) is even 
more solid, since tu cancels out in the ratio of U e 3 to U e 2- 



Conclusions 

We have presented a model of quark and lepton masses and mixings that is 
based on supersymmetric 5*0(10). The mass matrices arise from five Yukawa 
terms, each of which is of a type that arises in a straightforward way from 
simple tree-level diagrams. The model contains very few parameters, so few 
that there are nine predictions. 

Four of the predictions follow from the basic structure of the model and 
simple group-theoretic considerations. These are m u ~ (which corresponds 
to the observed fact that m u /m t ~ 10 -5 , much less than m^/mf, = 10~ 3 
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and m e /m T = 0.3 x 1CT 3 ); = mj (the well-known and successful £77(5) 
relation); m^/m® = |m°/m° (one of the Georgi-Jarlskog relations); and 
m Q d /m° s = Qm^/m ^ (the other Georgi-Jarlskog relation). 

Four other predictions are non-trivial quantitative successes of the model. 
They are a prediction of the charm quark mass that is correct within the 
present uncertainties; the prediction of very large mixing between and u T , 
consistent with the recent discovery of atmospheric neutrino oscillations; a 
prediction for V u b that slices neatly through the middle of the allowed region 
for the Wolfenstein parameters p and r/; and a prediction for sin 2 29 efl that 
is well inside the present allowed range for the small-angle MSW solution of 
the solar neutrino problem. 

The last prediction is of a presently unmeasured quantity: sin 2 29 eT = 
1.1 sin 2 26 efl . 

The model is eminently testable. The most dramatic test would be a 
measurement of 8 eT , which could come from improved proton accelerator 
bounds or perhaps from a muon collider. The next most important and 
fairly clean tests would come from a resolution of the solar neutrino problem 
and the consequent precise determination of 6 eil ; and from a measurement 
of the Wolfenstein parameters p and rj. Much more difficult, but hopefully 
achievable someday, would be precision tests of the predictions for m& and 
m c . This would require considerable knowledge of the spectrum of the spar- 
ticles, knowledge of tan (3 and the Higgs masses, and better determinations 
of a s , m t , and V cb Q. 
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Figure Captions 



Fig. 1: The diagram that generates the "e" entries of the quark and lepton 
mass matrices. 

Fig. 2: The diagram that generates the "<r" entries of the mass matrices L 
and D. 

Fig. 3: The diagrams that produce small masses for the first-family quarks 
and leptons. The family index j on the spinor field takes the values 2 or 
3, giving, respectively, the S and 8' terms of the mass matrices. Different 
SO (10) vector Higgs, denoted with superscript j, are exchanged in the two 
cases. 

Fig. 4 A comparison of the model with experiment for the Wolfenstein 
parameters (p, rf). The axes are p and r] multiplied by the central value 
of |si 2 Vrf,|. The central values allowed by the model lie on the bold dashed 
circular arc, cf. Eq. (6). The constraints following from \V u f,\, -B-mixing and 
e extractions from experimental data are shown in the lightly shaded regions. 
The experimentally allowed region is indicated by the heavily shaded central 
region. A typical unitarity triangle allowed by both data and the model is 
shown. 
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